IFUP-TH 76/94 



Running coupling constant and correlation length from Wilson 

loops 



Massimo Campostrini, Paolo Rossi, and Ettore Vicari 
Dipartimento di Fisica deU'Universita and I.N.F.N., 1-56126 Pisa, Italy 



Abstract 

We consider a definition of the QCD running coupling constant a{^) related 
to Wilson loops of size rxt with arbitrary fixed t/r. The schemes defined by 
these couplings are very close to the MS scheme (i.e. the one-loop perturbative 
correction to the coupling is small) for all values of t/r; in the t/r ^ oo limit, 
the "gg force" scheme is recovered, where the coupling constant is related to 
the quark-antiquark force. 

We discuss the possibility of applying finite-size scaling techniques to the 
Monte Carlo evaluation of a{fj,) up to very large momentum scales. We pro- 
pose a definition of correlation length, also related to Wilson loops, which 
should make such a computation feasible. 
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Typeset using REVTeX 



One of the most important goals of a non-perturbative approach to QCD such as the 
lattice is the determination of in units of a physical mass scale; this determination 
requires a very precise measurement of the running coupling constant in the MS scheme 
'^Ms(/^) high-momentum region, where perturbation theory is reliable and accurate. 

A determination of ajj^{fi^mzo) would be most welcome from a phenomenological point of 
view. Although considerable progress has been achieved in the last few years, especially in 
quenched QCD this task turns out to be quite hard, due to the difficulty of performing 
lattice calculations at large momentum scales. 

In the present letter we will address the two problems of: choice of an optimal for 
the lattice computation; determination of a(/i) up to the highest possible momentum scale. 

The requirements for an optimal choice of running coupling constant to be used in a 
lattice Monte Carlo simulation are: 

1) fast and accurate lattice measurement; 

2) proximity to a^js, to reduce systematic errors due to neglected orders in the pertur- 
bative conversion to the MS scheme; 

3) well-defined finite-size scaling properties. 

We shall now present a new definition of addressing the above points. 
A natural definition of running coupling is derived from the static quark-antiquark force 
F{r) by the relationship 

^W = -e.^, (1) 

where cp = {N'^ ~ 1)/ (2A^) = |. (We prefer to write agg with a momentum scale dependence, 
in analogy with aMs(yu), rather then with a length scale dependence, as often seen in the 
literature). In perturbation theory the qq and the MS couplings are related by a very 
small first-order coefficient , making the determination of ajj^ from agg by a perturbative 
redefinition of the coupling at large momentum very precise. 

The force F(r) between two static quarks separated by a distance r can be evaluated 
from rectangular Wilson loops W{r,t): 

V{r) = - \im - In W{r, t). (2) 

Recent developments in the calculation of agg were reported in Refs. 

On the lattice the implementation of this idea to extract (ygq{r) presents two major 
problems: 

1) it requires a t — > oo limit procedure in order to evaluate V{r); 

2) it requires control of the force from short to long distance in one simulation, with 
a <C r <C L holding for the whole range of physical r involved. This strongly limits the 
accessible values of r since the lattice size L/a available on today's supercomputers is not 
very large. 

In this letter we extract a{fi) from rectangular Wilson loops W{r,t) of size r x t with 
X = t/r fixed. This solves the problem (1) above, and opens the road to the use of finite 



2 



size techniques in order to reach very small distances without the necessity of very large 
lattices. Moreover, the close connection with the MS scheme is retained, as we shall see 
from a perturbative calculation. 
In weak coupling. 



X{r,t) 



92 \nW{r,t) 
drdt 



a 



CFA{x)—{l + 0{a)) , 



(3) 



where x = t/r and 



A{x) = - 



1 arctan(l/x) 
arctan(x) + h 



X 



x^ 



(4) 



The renormalization properties of the Wilson loop operator 
running coupling constant a^il/r) parametrized by a; = t/r: 



X{r,t) = -cfA{x) 



ax(l/r) 



allow the definition of a 



(5) 



Without loss of generality, we will choose r to be the shorter side of the loop, and therefore 
X > 1. Since 



lim = F{r), 

t — >oo 



(6) 



we have 



lim Q;^(l/r) = agg{l/r). 



(7) 



The perturbative expansion of xi"^, t) in dimensional regularization can be obtained from 
the corresponding expansion of the Wilson loop W{r,t) In the MS renormalization 

scheme 



X{r,t) = -CFA{x)^^i^ 



1 + 60 (ln(r/i) + R{x)) aMs(/^) + O (a^) 



(8) 



where 60 is the first coefficient of the /^-function: 

d 

1/22 4 \ 
bo = — ( —N --Nf), 
" 47r V 3 3 V 

and R{x) is a finite function of x: 



(9) 



R{x) = 7s - 1 + 



AttBo V 9 
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+ 



7rA{x) 



hix) + Ux) + 



N 



(I,(x) + h{x) + h{x)) 



(10) 
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The functions /i, I2, 13, h, h are given in the Appendix. They all vanish in the large-a; limit 
where the perturbative expansion of the static quark-antiquark force [0 is recovered. 

R{x) determines the first order connection between the couplings a^il/r) and ajj^{fi): 

«^(l/r) = a.(l/r) - bo R{x) aUl/r) + O (a^(l/r)) , 

where A^^ is the x-dependent A-parameter associated with the x-schemes 

In Table | we report R{x) for Nf = 0, 2, and 4, for various values of x. In quenched QCD 
{Nf = 0), R{x) turns out to be independent of A^; the values for Nf 7^ in the Table are 
computed for SU(3). As anticipated above, the first order perturbative coefficient relating 
to ajj^ remains small for all x > 1. 
On the lattice we can define the quantity 

, , , f r xr\ , ^ 

0(r,a;) = xc ' (12) 



a a 



where xc is the usual Creutz ratio 

(i,T\ 1 W{R,T)W{R+l,T+l) 
Xc{R, T) = - In ^(^^^+1)^(^+1^^)- (13) 

0(r, x) is a natural estimator of x('";2:r). By measuring 0(r, x) at different scales r keeping 
X fixed, one may easily extract the corresponding running coupling ^^.(l/r). 

We would like to conclude with some considerations concerning the possibility of em- 
ploying finite size techniques to reach very small distances, without the need of very large 
lattices. The power of this kind of techniques was recently illustrated by the exploration of 
the extremely small distance regime of two-dimensional spin models |[Tl[ . 



In the scaling region the following finite size scaling relations holds 

Uf3)^ML/^L)Um, (14) 

where ^ is a correlation length of the theory, and 

r, X, L) ~ /<^(x, L/r, L/^l) r, x, 00), (15) 

Notice that the existence of an additional scale r in 0(r, x) causes an extra dependence on 
the ratio L/r in the finite size scaling function f^. Finite size scaling functions like and 
may be reconstructed by performing simulations at relatively small lattices, by employing 
the scheme outlined in Ref. [^. In this context finite size methods should allow us to keep 
a <C r even when r is very small in physical unit. 

A crucial point in the realization of this program is the definition of a correlation length 
^ suitable for an accurate finite size scahng study. Once found such a scale, the study of 
the finite size scaling of should be relatively easy. Finite size studies of two-dimensional 



spin models [11,12| suggest that a suitable correlation length is one defined from the second 



moment of a correlation function. Such a correlation function should have appropriate 
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properties: it must renormalize multiplicatively with a renormalization factor independent 
of the distance, fall down exponentially at large distances, and of course be easily measurable 
in Monte Carlo simulations. We identified a possible candidate for pure gauge theory, which 
is easily constructed from Wilson loops. The correlation function 

Y(r t) = ^^'^'^ (16) 

is renormalized by a constant factor, since the divergence associated with the perimeter term 
of the Wilson loop operator cancels in the ratio; it is exponentially suppressed at large 
distances, due to the area law of pure gauge theory; it is constructed from Wilson loops, 
which are a standard "easy" observable of Monte Carlo simulations. 
From Y{r,t) we can define a second moment type correlation length: 

^2 l Jo°"drjrdtY{r,t)rt 

^ 2 J,°°drjrdtY{r,t) ' ^ 

where k is a free parameter (/t>l), which can be chosen to optimize the measurement. 
Assuming an exact area law for the Wilson loop, Eq. ([T?!) would give = l/cr, where a is 
the string tension. 

The measurement of ax{fJ^) in the above-mentioned scheme up to a large /i leads to a direct 
determination of the /3-function of the SU(3) lattice gauge theory and of the adimensional 
quantity ^Aj^. This quantity still needs to be converted to a more phenomeno logical scale, 
such as tq [|13|; but this is a rather minor problem, since it involves only a measurement of 



^/ro, which can be performed at the values of (3 of our choice. 

The program presented here is not specific to pure gauge theories, and both the definition 
of coupling constant and the finite-size scaling relationships hold unchanged for full QCD, 
with the inclusion of dynamical fermions. The definition of correlation length has to be 
changed; however a "physical" definition such as the inverse nucleon mass is perfectly viable 
in this case. 
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APPENDIX: 

hix) = ^Aix) - r| [fix) + fil/x)] - \nx] , 

^ , , 7rln2 ^ , , 1 2 r / X /. / M 

h{x) = — A{x) - 2^ ^ lai^) + 9i^/^)] ' 

h{x) = -r'^^[h{x) + h{l/x)], 
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h{x) = l{x) + ^-^, (Al) 
where x — t/r and 
f{x) — xarctan(a;) — - ln(l + x^), 

rx 

h{x) = f{x) (6 - 41nx) - (1 + x^)arctan^(x) + 2 dy [21ny(x - y) + f{x - y)] arctan(y), 

k{x) = -4/(x)/(l/x), 

/(a;) = — / / dzidz2 / c?p 

pZ2{.y'^ - + 2x2:i - 2y2:i) 
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TABLES 





TABLE L We report R(x) = InA^/Avr^ for Nf 


= 0, 2, 4 and for some 


values of x. 


X 


Nf = 


Nf = 2 


Nf = 4 


1 


-0.20758 


-0.27875 


-0.37270 


5/4 


-0.11177 


-0.18145 


-0.27343 


4/3 


-0.08975 


-0.15852 


-0.24930 


3/2 


-0.05567 


-0.12242 


-0.21053 








-D 18146 


7/4 


-0.02178 


-0.08554 


-0.16970 


2 


-0.00069 


-0.06189 


-0.14266 


3 


0.03187 


-0.02341 


-0.09637 


4 


0.03979 


-0.01317 


-0.08308 


10 


0.04557 


-0.00500 


-0.07175 


CX) 


0.04691 


-0.00324 


-0.06945 
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